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Abstract Social networks are discrete systems with a large amount of het- 
erogeneity among nodes (individuals). Measures of centrality aim at a quan- 
tification of nodes' importance for structure and function. Here we ask to 
which extent the most central nodes can be found by purely local search. We 
find that many networks have close-to-optimal searchability under eigenvec- 
tor centrality, outperforming searches for degree and betweenness. Searchabil- 
ity of the strongest spreaders in epidemic dynamics tends to be substantially 
larger for supercritical than for subcritical spreading. 
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1 Introduction 

The science of networks |36J is an interdisciplinary enterprise with application 
areas ranging from systems biology [3] and systems chemistry [39' to ecology 
[7] and the social sciences [3S]. Methods and models have been traditionally 
developed in the fields of graph theory [13j and algorithms [461IT2] . Dealing 
with large disordered and heterogeneous systems, statistical physics now has 
significant influence on network science [2,!15]. The analysis of networks fur- 
ther beneflts from the growing connections between statistical physics and 
computational complexity (35j . 

On the empirical side, advances in technology for observation and mea- 
surement of large systems and for compilation and storage of the resulting 
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data sets have triggered a revolution. In the social sciences, interaction net- 
works typically comprising less than 100 individuals used to be obtained by 
surveys [H] or by direct long-term observation of social groups, e.g. ^T\ . 
Nowadays digital communication by mobile phones 24, 4D] and e-mail [Ml 
140] and through web-based social media such as Facehook, Twitter and Flickr 
[22] automatically generates data sets covering millions of individuals. 

An important property of social networks is searchability, i.e. the pos- 
sibility of finding a target node by iteratively following connections (edges) 
and using information on properties of nodes in the direct neighborhood for 
deciding which connection to follow [Tll27|. This locally available informa- 
tion, such as the geographical location of an individual and other personal 
attributes, may be represented by placing nodes in a Euclidean space A 
searcher then iteratively chooses the edge that takes it closest to the target. 
The success and efficiency of such network navigation crucially depends on 
the distribution of long- and short-range connections in the network [?7] . 
Social experiments indicate that these conditions are fulfilled because per- 
sonal messages tend to reach the addressee in a few steps when iteratively 
forwarded by one person to a chosen acquaintance [551H71fH) . An informa- 
tion theoretic approach characterizes the reachability between nodes by the 
minimal description length of the search path [^H5] . 

A related task of large practical relevance is the search for the most 
important nodes in a given context. These are, for instance, those that would 
contribute most to the spreading of an epidemic contagion and thus should 
be vaccinated first [4T | I11 ) . In spreading and other dynamical processes, the 
importance of a node strongly depends on the network context |261I30[ . It may 
be captured by centrality measures such as the number of neighbours, the 
number of shortest communication paths running through, or the number 
of walks emanating from the node [48ll31j. Given a centrality measure, a 
landscape is obtained by labeling each node with its centrality value [SUlli]. 

Here we investigate which network structures and centrality measures give 
rise to searchable landscapes where the most central nodes are eventually 
reached by iteratively jumping to the most central neighbour, i.e. by local 
search. We quantitatively answer this question by introducing the smoothness 
as the expectation value of the centrality eventually encountered, normalized 
by the maximum centrality in the whole network. 

In practical scenarios, local search can be applied efficiently only when 
the centrality of a node can be obtained locally at that node. Eigenvector 
centrality, for instance, is of global nature. The explicit and exact computa- 
tion of eigenvector centrality necessarily involves the whole network and then 
yields the centrality values of all nodes. Once all these values are known, the 
maximum is readily identified. In this case there is no need for local search 
on the centrality landscape. 

Local search is useful when the network as a whole is not known but 
centrality of an encountered node can be obtained indirectly, e.g. by local 
measurement of system dynamics. The degree of a node is readily obtained 
by measurement of the local density from diffusion. The stationary fraction 
of random walkers at a node is proportional to the node's number of neigh- 
bours. Eigenvector centrality is measured approximately as the frequency 
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with which a node is involved in critical spreading or percolation clusters 
[30] , We term this node's property the spreading centrality. The smooth- 
ness of the corresponding landscape quantifies the success of identifying the 
globally strongest spreaders by iteratively building a path towards stronger 
spreaders. 



2 Definitions and Methods 

2.1 Networks and landscapes 

A network (also called a graph) is a tuple G — {V, E) with V a finite set and E 
a set of unordered tuples in V. The elements of V are called nodes, elements of 
E are called edges. For a node v €V oi a. network (V, E), the neighbourhood 
N{v) is the set {w € V : {v, w} € E}. The closed neighbourhood of v is 
N{v) = N{x) U {v). 

If G = (y, E) is a network and / : — K is an arbitrary mapping, 
then L — {V,E,f) is called a landscape (over G). The max-neighbourhood 
of a node is the set Nmax{v) — {w N{v) : f{w) = maXxgAr(t,) f{x)}. A 
node V G V is [strict] local maximum of the landscape if, for all w £ N{v), 
f{v) > f{w) [f{v) > f{w)]. 



2.2 Search dynamics and smoothness 

An adaptive walk is a local search dynamics. At each time step t, a transition 
from the current node i; to a neighbour is proposed by drawing w G N{v) 
uniformly. If f{w) > /(w), w is accepted as the node for the next time step, 
otherwise the search stays at v. Formally, an adaptive walk on G is a homo- 
geneous (time- independent) Markov chain with state set V and transition 
probability 

^ I .(W if -e^W-d /(.«)>/(.) 
^ lO otherwise ^ ^ 

for w V, and Tr{v — > v) given by normalization of probabilities. An adaptive 
walk is equivalent to kinetics at zero temperature where — / plays the role 
of energy and the network encodes the allowed transitions between configu- 
rations (nodes). If w is a strict local maximum of L, then 7r(u v) — 1, so v 
is an absorbing state for the adaptive walk. 

A gradient walk is similar to an adaptive walk, but concentrating on /- 
maximal neighbours. The dynamics proceeds to a neighbour with the max- 
imal value for / in the neighbourhood, provided that this maximal value is 
strictly larger than the one at the current node. A gradient walk is thus a 
Markov chain with transition probability 

^w) = i ^ and fiw) > f{v) 

I otherwise 
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for w ^ V, and n(v — > v) given by normalization. A node v is an absorbing 
state of the gradient walk (7r(u v) = 1), if and only if u is a local maximum 
of the landscape. 

Now we characterize the landscape by the success of the search dynamics. 
By (/) {t) we denote the expectation value of / under a given search dynamics 
with the uniform distribution on the node set V as initial condition. For finite 
time i, we measure the t-smoothness of L, 

St(L) = (/>(t)/max/(«) (3) 

to quantify how close the search dynamics approaches nodes with maximum 
centrality. The long-term success of the search is quantified by 

s{L) = lim st{L) , (4) 

>C30 

the smoothness of L. The limit in Eq. Q exists because none of the search 
steps decreases the /-value (so this holds also for the expectation value) and 
/ is of finite support, thus upper-bounded. For non-negative functions / — 
such as the node centralities considered in the following — the smoothness 
takes values in the unit interval. In particular, s{L) = 1 indicates that the 
search is certain to reach a node with maximum /-value, regardless of initial 
condition. 

The opposite case, s{L) < 1, occurs exactly when the Markov chain has 
an ergodic set including a node on which / is not maximal. For the gradient 
walk, this means that there is a local but not global maximum of / because 
here the ergodic sets are exactly the singletons formed by local maxima. In 
the adaptive walk, each strict local maximum v* gives rise to an ergodic set 
{v*}; further ergodic sets are formed by interconnected nodes of the same 
centrality value, all without neighbours of larger centrality. 



2.3 Centrality measures 

The degree centrality of a node v is defined as the number of neighbours 

div) = \Niv)\ . (5) 

The degree is a purely local measure of node importance, taking into account 
only the neighbouring nodes of the node considered. The degree distribution 
P Tn U {0} ^ M is defined by 

p(« = . (6) 

The fc-core of a network G is the largest subnetwork of G in which all 
nodes have degree at least k. It is computed by iteratively removing nodes 
with less than k neighbours until no such nodes remain in the network. The 
removal of a node deletes all its edges so the loss of node v may suppress the 
degree of the neighbours of v below k as well. Therefore the node removal is 
done iteratively. 
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The shell index h(v) of node v is the largest number k such that v is in 
the fc-core of the network. A large shell index indicates that the node is part 
of a densely connected subnetwork. 

The eigenvector centrality is obtained by finding a non-negative centrality 
value e(v) for each v ^ V to solve the set of equations 



with Amax G R the maximum value for which such a solution exists. The 
underlying idea is that the importance of a node is obtained self-consistently 
as the sum of its neighbours' importances discounted by a node-independent 
factor A~3jx. Equation Q can be written as the eigenvector equation Amaxe = 
Ae for matrix A, hence the name. The adjacency matrix A has entry Uvu, = 
1 when {v, w} is an edge, and = otherwise. The Perron-Frobenius 
theorem guarantees that the solution is unique up to a common scaling factor 
if the network is connected, i.e. for any two nodes v and w, there is a path 
between v and w. 

The betweenness centrality describes node importance as the property of 
being contained in shortest paths between other nodes. We denote by auw 
the number of shortest paths between nodes u and w, by cr„u,(w) the number 
of such paths passing through node v. Then 



is the betwenness centrality of node v. Computation of betweenness centrality 
for all nodes in an unweighted network {V,E) requires 0(1^1 computa- 
tional time steps and OdV^I + \E\) memory. 

2.4 Spreading centrality 

Bond percolation |23) is a theory for the description of spreading of material 
and information in disordered media. On arbitrary networks, the statistics 
of clusters generated by bond percolation coincides with those of the basic 
Susceptible-Infected-Removed (SIR) model of epidemic spreading [38ll25j. 

A realization of bond percolation on a network G — (V, E) at parame- 
ter /3 is a network R = {V, E') obtained as follows. For each edge e & E 
independently, include e £ E' with probability /?, omit e with probability 
1 — /3. Such a realization defines a partition of the node set V into clusters 
Ci, C2, . . . . Each node is in exactly one cluster. Two nodes v and w are in 
the same cluster if and only if E' contains edges by which one can walk from 
V to w. By zp{v) we denote the expected size of the cluster containing v in 
bond percolation at parameter /3. The spreading centrality (t>p{v) is obtained 
by normalization according to 




(7) 



weN(v) 





(9) 



max^gy zi3{w) - Zfi 



6 



with the mean value = Ylw^v simulations, Zj^^v) is obtained 

from 10^ independent realizations of bond percolation for each network and 
value of /3. 

For each node v, monotonically increases with /3. However, the 

increase may differ from node to node. Therefore the ranking of nodes with 
respect to spreading centrality substantially changes with varying f3 [3D]. 
Then the landscape of spreading centrality L — {y,E,(j)p) is qualitatively 
different for different values of (3. 



2.5 Random networks 

Classic random models of networks are by Gilbert [20 and Erdos and Renyi 
[19j . Gilbert's random graph model has two parameters, being the num- 
ber of nodes n S N and the edge probability p e [0, 1]. Among all graphs 
on n vertices, each graph with rj edges is obtained with probability p^{\ — 
p^jn{n-i)/2-r} ^ In Other words, a realization of the model decides independently 
for each of the n{n — l)/2 possible edges if the edge is present (with proba- 
bility p) or absent (with probability 1 ~ p). Throughout this contribution we 
refer to Gilbert's model as random graph. 

The model by Erdos and Renyi is defined as the uniform distribution on 
all graphs with n nodes and exactly iV edges, n e N and G {0, n{n — 1) /2} 
being the model parameters. Asymptotically, for n — > oo, the models by 
Gilbert and by Erdos and Renyi become arbitrarily similar when taking rj = 
pn{n — l)/2. 

Another, more versatile form of statistical ensembles are those with a 
given degree sequence. The degree of each of the n nodes is a parameter 
value of the ensemble. The ensemble consists of all graphs with nodes hav- 
ing the prescribed degrees, having uniform probability. Here we use Markov 
chain Monte Carlo sampling with edge switching |43j in order to randomize 
networks under conservation of degrees. In a step of edge switching, two node- 
disjoint randomly chosen edges {v,w} and {x,y} are drawn (proposal) and 
replaced by edges {v, x} and {w, y} (acceptance) unless these edges already 
exist (rejection). In order to obtain a randomization of a network with n 
nodes, we run the Markov chain until switching has been performed times 
(counting only accepted steps) . The configuration model [M] is an alternative 
method for the same purpose. 



2.6 Stochastic network growth 

Let us introduce four procedures for building up networks by iterative ad- 
dition of nodes and edges. The following scenario is common to all these 
procedures. The network is initialized with m nodes fully interconnected, i.e. 
all m{m— l)/2 edges are present, where to is a parameter. Then in each step 
t of growth, a subset S of size m is chosen stochastically from the current 
node set; a new node v{t) and to new edges are added to the network, where 
each edge connects one node in S with v{t). The growth step is iterated until 
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the network reaches the desired size. Each specific procedure is defined by 
the stochastic choice of the set of nodes that the new edges are attached to. 

The preferential attachment rule [51 l32l[TC] (also termed cumulative advan- 
tage im) adds edges preferably to nodes having large degree already. Starting 
with an initially empty set S, a node v is drawn from the distribution 

/ N d{v) + a 

= , , r (10) 

E^ev\sid{w) + a) 

on V \ S. Then this node is included in S. This drawing of nodes without 
replacement is repeated until \S\ = m. The parameter a is a constant bias 
with m < a < oo. Uniform attachment, called model A in refs. OH], is 
obtained in the limit a — > oo. Here S{t) is drawn uniformly from the set of 
all m-node subsets. 

The edge attachment rule by Dorogovtsev et al. 'il7] is defined specifically 
for m = 2 edges to be added per node. Rather than selecting nodes separately, 
the set S is determined by drawing an edge e G E uniformly and setting 
S = e. Thus each new node is attached to the end nodes of an existing edge, 
thereby forming a triangle. 

In the deactivation model by Klemm and Egufluz ^291I281[T5| , the set of S 
of nodes receiving edges — thus called active — evolves over time. Upon initial- 
ization, S comprises all m nodes present initially. In each growth each step, 
after adding a node v and its m edges, v is included in S. From the |5| = m+l 
then contained in S, one node u is chosen for removal (deactivation) by the 
distribution 

/ X d(u)~^ . . 

E^es^W 1 

The expected degree distribution from the deactivation model and from 
edge attachment decay as a power law P{k) ~ k~'^ with 7^3. Preferential 
attachment yields P{k) ^ k^^ with 7 = 8 + a/m. The asymptotic degree 
distribution from uniform attachment decays geometrically. The deactivation 
model generates networks with average distance between nodes increasing 
linearly with size; the other three rules generate networks with distances 
increasing logarithmically or sublogarithmically. 

Note that all these procedures generate networks with shell index h{v) = 
m for all nodes v, as shown by induction over the number of nodes. The 
initial condition is a fully interconnected network of m nodes. After adding 
a node and m links, each node v now has exactly d{v) — m neighbours and 
therefore h{v) = m. After an arbitrary number of node additions, each node 
has at least m neighbours so that the m-core is the whole network. In the 
computation of the m + 1-core, nodes with degree m or smaller are iteratively 
removed, which amounts to disintegrating the whole network in node order 
similar to its generation. Thus the m + 1-core is empty. Each node has the 
same shell index m, which is the global maximum. It follows trivially that 
all these networks have smoothness 1 with respect to shell index. 
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Fig. 1 Smoothness of centrality landscapes on stochastically grown networks. 
Symbols distinguish centrality measures (see legend). Large open symbols are for 
search by adaptive walk, small filled symbols for gradient walk. Panels (a-f) dis- 
tinguish procedures for network generation, (a) preferential attachment with bias 
a — 0, degree exponent 7 = 3 (b) preferential attachment with bias a = — m/2, 
7 = 2.5. (c) edge attachment, (d) uniform attachment, (e) deactivation model, (f) 
random graph with edge probability p = 4/(n— 1). Error bars indicate standard 
deviation over 10 independent realizations of network generation. In the growing 
networks (a)-(e), m = 2 edges per node are attached in the growth process; in (f) 
edge probability is p = 4/(n — 1), resulting in asymptotic average degree of 4 in all 
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Table 1 Smoothness of centrality landscapes for six empirical networks. In each 
column, the first value gives the smoothness for search by adaptive walk. The 
second value is for search by gradient walk. For each network, the first (upper) row 
of smoothness values is for the original network. The second ro w is f or randomized 
surrogate networks with the same degree sequence (see section 2.5 1; each value is 
the mean of the smoothness values of 10 independent randomizations . Networks 
are (a) the e-mail cont acts from University Rovira i Virgili, restricted to the largest 
connecte d co mponent [24j : (b) jazz bands connected by an edge if they share a 
musician [21]; (c) authors in the cond-mat e-p rint archive (arXiv) where {«, w} is an 
edge if v and w have co-authored a paper [37] , updated version including data until 
March 2005; (d) users of the Pretty-Good-Privacy algorithm for secure information 
interchange ;8j (e) Internet at Autnomous Systems level, snapshot taken on July 
06, 2006. (f) electric power grid with generators, transformers and substations as 
nodes, edges being high- voltage transmission lines US]. Network data (a), (b) and 
(d) were downloade d from h ttp : //deim.urv . cat/~aarenas/data/welc ome .htm, ; 
(c), (e) and (f) from jhttp : //www-personal .umich. edu/~mejn/netdata/| 



network 




degree 


eigenvector 


shell index 


betw.ness 


(a) e-mail 


1133, 5451 


0.882 
0.859 


0.949 
0.892 


1.000 
0.905 


1.000 
0.939 


0.999 
0.992 


1.000 
0.992 


0.867 
0.791 


0.908 
0.851 


(b) jazz 


198, 2742 


1.000 
1.000 


1.000 
1.000 


1.000 
1.000 


1.000 
1.000 


1.000 
1.000 


1.000 
1.000 


0.830 
1.000 


0.934 
1.000 


(c) cond-mat 


40421, 175692 


0.766 
0.951 


0.773 
0.967 


0.921 
0.957 


0.921 
0.971 


0.791 
0.995 


0.794 
0.995 


0.627 
0.926 


0.671 
0.955 


(d) pgp 


10680, 24340 


0.569 
0.832 


0.607 
0.853 


1.000 
0.861 


1.000 
0.870 


0.848 
0.898 


0.852 
0.898 


0.736 
0.813 


0.783 
0.846 


(e) Internet 


22963, 48436 


0.973 
0.869 


0.981 
0.872 


1.000 
0.892 


1.000 
0.892 


1.000 
0.895 


1.000 
0.895 


0.996 
0.888 


0.997 
0.889 


(f) power grid 


4941, 6594 


0.401 
0.474 


0.410 
0.495 


0.230 
0.623 


0.233 
0.630 


0.651 
0.970 


0.651 
0.970 


0.411 
0.395 


0.437 
0.426 



3 Numerical results 



We first study smoothness of centrality landscapes for stochastically gener- 
ated network instances, as described in sections |2.5| and [2.6| The dependence 
of smoothness on the number of nodes is displayed in Figure [T] Smoothness 
depends both on the type of network and the centrality measure under con- 
sideration. Except for the random graphs and the networks grown by the de- 
activation model, eigenvector landscapes have a maximum smoothness s = 1 
with deviations only for small networks. Betweenness centrality landscapes 
have smoothness close to 1 in networks grown by edge attachment. For all 
other combinations of network generation and centrality measure, a decrease 
of smoothness with the size of the network is observed. In those cases, cen- 
trality landscapes become less searchable with increasing size. 

The aforementioned results are qualitatively the same for both types of 
local search. Except for the networks from uniform attachment, smoothness 
values from adaptive and from gradient walks are also quantitatively the 
same up to statistically insignificant deviations. 
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Fig. 2 Smoothness of spreading centrality landscapes under adaptive walks (open 
squares) and gradient walks (filled squares). Vertical dashed lines indicate the pa- 
rameter value /3c =^1/Amax as an estimate of the transition point from local to 
global spreading [10]. Curves with number symbols t £ {1,2,3} show S(L), the 
t-smoothness of the adaptive walk at time t, cf. Equation ([3|. See the caption of 
Table [l] for details on the networks. 



Now we turn to centrality landscapes based on empirical networks. Table[T] 
provides an overview of smoothness values for four social networks (a-d) and 
two technological networks (e-f) further described in the caption. We find 
that eigenvector centrality induces the smoothest landscapes on the social 
networks and on the Internet snapshot (e). On the power grid (f), however, 
eigenvector centrality reaches lower smoothness than the other three cen- 
tralities. This network is also the minimum of each of the eight smoothness 
columns, i.e. the power grid yields the lowest smoothness among the networks 
under a given centrality measure and walk type. 

Randomization of these networks under conservation of the degree se- 
quence drastically changes the smoothness values. The searchability of cen- 
trality landscapes depends on properties of the empirical networks beyond 
the degree distribution. 

For landscapes of spreading centrality (see section 2.4), the smoothness 
as a function of the spreading parameter /3 is plotted in Figure [2] For small 
/3, smoothness values are close to those of the corresponding degree land- 
scapes, see Table [l] This finding is in accordance with the large rank order 
correlation between spreading centrality and degree below the transition |30j . 
In the transition regime (/3 w /3c), smoothness increases with /3. The rank 
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order correlation between degree and spreading centrality is also large in the 
supercritical regime. Here, however, smoothness for spreading centrality does 
not coincide with smoothness for degree. The /3-dependence of smoothness 
varies across networks in the supercritical regime. 

Above the transition, the i-smoothness values st{L) for small number 
of time steps t strongly increases, indicating an accelerated success of the 
search for more central nodes. Typically in the supercritical regime, t = 2 
or t = 3 steps are sufficient for the adaptive walk to saturate, i.e. to reach 
St{L) « s(L). 



4 Estimates for a given degree distribution 

Let us estimate the smoothness of degree centrality landscapes on the basis 
of the degree distribution alone. To this end, we perform an annealed ap- 
proximation. At each time step, the network is drawn uniformly at random 
from the set of networks with the given degree distribution. 

Prom the given degree distribution P{k) we obtain the probability of 
encountering a node of degree k by following a uniformly chosen edge as 

P'{k) = kP{k)/k (12) 

with the average degree k = Y1^=q kP{k). We use the cumulative of P' being 

oo 

Q'{k) = Y,P'{k). (13) 

i—k 

We define x{k,t) as the fraction of searches being active at a node of 
degree k at time t. A search is active at initialization and turns inactive 
when reaching a node that is a strict local maximum w.r.t. degree. By y{x, t) 
we denote the fraction of searches being inactive at a node of degree k at 
time t. We have Yl'k'=o^(^^^) + y{k,t) = 1 as a normalization at all times 
t. Initially, all searches are active and distributed uniformly across nodes, so 
x{k, 0) = P{k) and y{k, 0) = for all k. 

If all neighbours of a node with degree k have degree strictly smaller 
than k, the adaptive walk stays there, so the search becomes inactive. This 
happens with probability 

a{k) = [l-Q\kt . (14) 

Thus we have 

y{k,t+l)=y{k,t)+a{k)x{k,t) (15) 

Otherwise, the searcher eventually finds a neighbouring node with degree 
at least k and jumps there. In order to simplify the equations, we assume 
that this jump happens immediately in one time step. Compared to the real 
dynamics, where one or several rejections may occur before the jump, this 
simplification modifies the transients but not the asymptotic distribution of 
searches across degrees. When following a uniformly chosen edge, we find a 
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Fig. 3 Smoothness estimates s* compared to the real smoothness for adaptive 
walks on degree centrality landscapes. Symbols distinguish networks as follows. 
Squares: 10 independent realizations of random graphs with n — 1000 nodes and 
parameter p — 4/(n— 1). Cycles: Same as for squares but n = 9000. Stars: empi rical 
networks randomized by switching under conserved degrees, cf. subsection |2.5[ The 
abscissa of each star is the smoothness averaged over 10 independent realizations of 
randomization (see also Table [T]). Error bars indicate the standard deviation over 
the 10 realizations. 



node of degree k with probability P'{k)/Q'{l) under the constraint k > I. 
Therefore the active searches propagate as 

k 

x{k,t + l)=Y,i^-amP'ik)/Q'il)]xil,t) . (16) 



When iterating Equations (15 15) over time steps, the fraction of active 
searches converges to zero. The estimate s*{L) of the smoothness is obtained 
from the asymptotic distribution of inactive searches 

CO 

s*(L) = hm Vy(fc,i)A:/fc^ax (17) 

5-00 ^ — ^ 
k=0 

with the maximum degree /cmax- 

Figure [3] provides a comparison between actual smoothness value s and 
the estimates s* for several networks. For realizations of random graphs, 
the degree distribution provides almost complete information on smooth- 
ness. The typical difference in smoothness s across realizations is significantly 
larger than the deviation of the estimate s* from the true value. Applied to 
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randomizations of empirical networks, there is reasonable agreement between 
s and s* in most of the cases. 



5 Discussion and outlook 

We have defined a framework for assessing the searchability of centrality 
landscapes arising from a network. The smoothness of a network combined 
with a centrality quantifies the extent to which nodes of large centrality are 
found by searching the network locally. Local search means that the sequence 
of nodes visited is a walk following edges on the network. 

This local perspective is motivated by the limited information on the 
network. Rather than knowing the whole network at the outset, a searching 
agent explores the system by iteratively following connections. In the present 
framework, we have assumed that the centrality value of each node becomes 
available ad hoc when encountering the node. This is the case when the 
centrality values themselves involve local information only, e.g. the degree. 
Alternatively, the centrality of interest results from local measurement on 
dynamics taking place on the network. We have investigated the spreading 
centrality as an example. 

Analytic insight into smoothness s of centrality landscapes is desirable, 
e.g. lower bounds on s depending on network properties. As a step in this 
direction, we have defined rate equations for estimating smoothness of degree 
landscapes in an annealed approximation. 

Numerically we find that eigenvector centrality typically generates maxi- 
mally (s = 1) or almost maximally smooth landscapes. As a heuristic expla- 
nation, smoothness arises due to the summation in the eigenvector equations 
that amounts to an averaging over neighbouring nodes. On the stochasti- 
cally generated networks lacking small- world property |28j and the extremely 
sparse electric power grid, none of the centrality measures under considera- 
tion generates smooth landscapes. Small neighbourhoods and long distances 
in search space tend to create obstacles to local search. 
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